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A. Taylor 

Abstract. Our main purpose is to give multiple examples for using the 
available implementations for computing the normalization of an afflne 
ring, computing the minimial generators of the normalization as an al- 
gebra over the original ring and integral closures of ideals. Some such 
examples have been published for Singular, but not for Macaulay 2 
and we present both in this paper. We also briefly describe the imple- 
mentations. 



1. Introduction 

In this note we describe the available methods for computing the nor- 
malization of affine rings, that is reduced, finite, fc-algebras, where k is a 
computable field, and the implications for computing the integral closure of 
ideals. Our goal is to keep it simple. Much more is known for computing 
the normalization with assumptions on the ring or the ideal, as well as for 
the normalization of discrete valuation rings. We only address what can be 
done in the most general setting of affine rings. 

Our main goal is to give examples so that users can fully exploit the im- 
plementations. We also briefly describe which algorithms have been imple- 
mented, some implementation issues and how to use these implementations 
to compute the integral closure of small ideals. All computations listed in 
this paper were done on a Pentium III, 600 MHZ machine with 256 MB of 
RAM. 

At the time of this writing there are two implementations of normal- 
ization. One implementation is in Macaulay 2 and the other is in 
Singular |10j . Both implementations are based on the information given 
in de Jong's paper |12j . The philosophical approach of which is seen in Vas- 
concelos' papers |16j and |18j and a more detailed analysis of this approach 
is given in [2]. A brief description as well as other examples for Singular 
are given in W. Finally, there will be chapter on normalization, comput- 
ing it in Singular, and related issues in the upcoming book of Greuel 
and Pfister jH]- All of what is published is based on the implementation 
in Singular, where as this is the first document detailing what is done in 
Macaulay 2 as well. 
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Let A = R/ 1 he an affine ring. In the following we use the language 
of Matsumura 1131 . If A is a domain and Q its ring of fractions then the 
integral closure A of A is the set of all elements in Q that are integral 
over A and A is integrally closed if A = A. An affine ring A is normal 
if Ap is an integrally closed domain at every prime p E Spec(^). If A 
is normal it is well known that A = A/ Pi x ■ ■ ■ A/Pt where Pi, ... i-^ are 
the minimal primes for A and A /Pi is an integrally closed domain for 1 < 
i < t. Since the implementations deal with affine rings in general this 
decomposition is often referred to as the normalization of A rather than 
the integral closure. Macaulay 2 refers to the computation as integral 
closure where as Singular and all of the publications associated to its 
implementation call the computation normalization. To emphasize that we 
do not require the ring to be a domain, we also refer to the computation as 
normalization. 

Serre's criterion for normality implies that an affine ring A is normal if 
and only if it is both S2 and Ri |131 Theorem 23.8]. An affine ring is ^2 
if Ap is Cohen-Macaulay for p G Spec{A) and codim(p) < 2, and has depth 
greater than or equal to 2 for p £ Spec(^) and codim(p) > 2. An affine ring 
is Ri if it is regular in codimension one. 

The philosophy of the approach is to enlarge A recursively inside of Q 
until the normalization of A is obtained. The extension Vasconcelos uses 
in |16| . jl8j is Hom(Hom(L, A), Hom(L, A)), where L is the Jacobian ideal 
of A. The following key theorem in de Jong's paper |12j is originally due 
to Grauert and Remmert [HI PP- 220-221], ,7, pp. 125-127] and describes 
the ring de Jong uses for this extension. The non-normal locus of A is the 
set NNL C Spec(A) such that for p G NNL, Ap is not an integrally closed 
domain. 

Theorem 1.1. Let A = R/I he an affine ring and J an ideal of A. 
Assume that the ideal J contains a non-zero divisor, and has the following 
property: 

NNL C T/(J), 

where V{J) = {p € Spec(74)| J C p} denotes, as usual, the zero set of J. 
Suppose moreover that J has the property 

(1.2) HomA( J, J) = HomA( J, A) D A. 

Then one has the following normality criterion: 

A = IIomyi(J, J) if and only if A is normal. 

The implementations in both MACAULAY 2 and Singular are split into 
two phases. First they find an ideal J such that the non- normal locus of 
A is contained in V{J) and satisfying Equations (|1.2|) . In the second phase 
IIom( J, J) is computed as a ring. 

For the first phase, Theorem 2.2 in de Jong's paper |12J establishes 
that any radical ideal containing a non-zero divisor satisfies Equation (|1.2|) . 
Also, it is well known that if the Jacobian ideal L of ^ is not contained in 
Q S Spec(A) then Aq is regular, so the non-normal locus of A is contained in 
V{L) = Hence \/L satisfies Theorem ll.il Computing L can be very 
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complex in terms of space and time, so sometimes the programs may use 
VT for f £ L non-zero which also satisfies the hypotheses of Theorem 11.11 

For the second phase, there are two facts that are very helpful. The first, 
is the well known fact that if / is a non-zero divisor in A, then (/J : J) = 
/Hom(J, J). The second, is the presentation of Hom(J, J) as a ring that is 
due to Catanese ^ and is given in de Jong's paper [12] Theorem 3.1]. Let 
vi, . . . ,Vn he a module basis for Hom( J, J). Map S = k[Xi, . . . , Xn] onto 
Hom( J, J) by sending Xi to Vi for 1 < i < n. Then as a ring Hom( J, J) is 
defined by two sets of relations. There are linear equations which are the 
syzygies for the module, ^ aiXi = and quadratic equations which come 
from the fact that Hom(J, J) is a ring so there exists Pijk, 1 < k < n such 
that ViVj = ^PijkVk- It is a theorem of Catanese stated clearly in l2! and 
|12j that the ideal generated by these two sets of relations is the full kernel 
of the map, thus inducing an isomorphism. 

The implementations in Macaulay 2 and Singular were completed 
independently, but the pseudo-code for both is very close to the presentation 
in [2j so we will not include it here. 

The casual user looking at the code in either program will see relatively 
complicated programs which is due to the recursive nature of the algorithm. 
The program must track the proper extensions and in the case of reduced 
rings, it must track the splitting of the ring into domains. Often for an 
ideal J as in the theorem, while Hom( J, J) is a proper extension, it is not 
the normalization yet so the process must be repeated many times. It is 
unknown, given the ring, the number of extensions that must be computed. 
Moreover, experimental experience suggests that there are choices available 
that could reduce the number of extensions that need to be computed, but 
knowing when to use them is not always a priori clear. There is a new 
algorithm of Vasconcelos' for which the number of recursions is bounded, 
but this algorithm is not yet implemented in any of the systems |17j . 

2. Examples 

For the examples we will focus on a hypersurface given to us by Craig 
Huneke, a union of straight lines given in the Singular example files and 
an example of Huneke's given in each of |2], |15) . |16j . and |18) . Also the 
first and third examples were chosen because one satisfies Serre's condition 
S2 but does not satisfy Serre's condition Ri and the other is the opposite. 

For each of the examples we give the input and output of both programs 
so that the user can simply type in the same information and see the same 
results. Occasionally, due to the width of the text we have had to slightly 
alter the output for it to fit in the space allotted. We hope that this is also 
enough information for the user to then perform similar computations on 
their own examples. Also, when possible we include the computation times 
in CPU seconds. If no time is listed it is because the time is negligible. 
We include the times because they are indicative of what can happen on 
larger examples and because they give more meaning to the times listed 
after Example 12.61 However, when taken by themselves the times may not 
have much meaning since even on the same machine times can vary up to a 
few seconds. 
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Example 2.1. Let A = <Q[x,y, z]/{x^ — — 'ip'z^). This ring is not 
Rl, but it is S2. Therefore, computing the normalization in the example 
forms the desingularization of A in codimension one. The computation in 
Macaulay 2 follows. 

11 : A = QQ[x,y,z] /ideal (x"6-z"6-y"2*z~4) 
ol = A 

ol : QuotientRing 

12 : time IntegralClosure (A, Variable => symbol a) 

— used 1.67 seconds 



QQ [a , a , X, y, z] 
7 6 

o2 = 

2 2 2 2 2 

(x -az, ax-az, a -ax, a -y -z) 
6 6 7 6 7 7 

o2 : QuotientRing 

Com,puting the normalization of an affine ring requires new variables unless 
the ring is already normal. Macaulay 2 gives the user the opportunity to 
specify which letter to use for this new variable. This is done via 

Variable => symbol a. 

// this is not included, then Macaulay 2 uses w by default. Using a letter 
to create the normalization which is a variable in the definition of R can 
lead to problems with using the normalization, so we recommend supplying 
a variable not used to define R. 

For Singular, we must first load the normalization library. Singular 
then lists the libraries loaded and we use elipses to indicate that more is 
printed on that line, but we have removed it to keep from having to make 
multiple lines due to the size of the paper. Also, we use nor to name the 
output of normal(I) to match the comments that singular outputs as part of 
this function. We then follow the other commands they list for continuity. 

> LIB "normal. lib"; 

// ** loaded /usr/local/Singular/2-0-3/LIB/normal.lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/lmoether.lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/primitiv.lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/random.lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/matrix.lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/ring. lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/inout . lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/presolve . lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/elim. lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/primdec . lib . . . 
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// ** loaded /usr/local/Singular/2-0-3/LIB/sing. lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/poly . lib ... 
// ** loaded /usr/local/Singular/2-0-3/LIB/general . lib ... 

> ring R = 0, (x,y,z) ,dp; 

> ideal 1= x6-z6-y2z4; 

> list nor = normal (I); 

// 'normal' created a list of 1 ring(s) . 
// nor [1+1] is the delta-invariant in case of choose=wd. 
//To see the rings, type (if the name of your list is nor) : 
show( nor) ; 

//To access the 1-st ring and map (similar for the others), 

type: def R = nor[l]; setring R; norid; normap; 
// R/norid is the 1-st ring of the normalization and 
// normap the map from the original basering to R/norid 

> timer=l; 

//used time: 0.82 sec 

> show (nor) ; 

// list, 1 element(s): 
[1] : 

// ring: (0) , (T(l) ,T(2) ,T(3) ,T(4) ,T(5) ) , 

(a(l,l,l,l,l),dp(5),C); 

// minpoly=0 
// objects belonging to this ring: 
// normap [0] ideal, 3 generator (s) 

// norid [0] ideal, 4 generator (s) 

> setring S; norid; 
norid [l]=T(4)-2-T(l)*T(5) 
norid[2]=T(l)*T(4)-T(3)*T(5) 
norid [3] =T(2) -2+T(3) "2-T(5) "2 
norid [4] =T(1) "2-T(3) *T(4) 

The the output of the command show(nor); states that the normaliza- 
tion of R/I requires only one ring because R/I was a domain. If we de- 
note the normalization of R/I by S/L then the next line tells us that S = 
Q[T(1), T(2), T(3), T(4), T(5)], each variable has degree 1 and the ordering 
is degree reverse lexicographic. The next two lines tell us that in Singular 
there are two other objects assigned to this ring. The first defines the nat- 
ural map from R/I into S/L and the second is the the defining ideal L of 
the integral closure. To see generators of L we execute the two commands 
setring S; and norid;. 

The information stored in the outputs of the two programs is fundamen- 
tally the same, but the presentation of the output in the two programs is 
a little different. Macaulay 2 preserves the fact that A C Hom( J, J) and 
thus computes a presentation that includes those variables from A which 
still contribute to a minimal presentation. In contrast, Singular uses all 
new variables and in this example T(l) = x,T(2) = y,T(3) = z,T(4) = 
aQ,T{5) = aj. The fractions computed in Example 12.41 can be used to find 
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this correspondence. Macaulay 2 allows the user to choose the name of 
the new variables, but Singular does not have such an option. When Sin- 
gular computes normal{I) it gives information about the number of rings 
in the decomposition and stores this information in the list we called nor 
above and then a few commands are used to display this information. In this 
example there is one ring, but in the next example when the input is not a 
domain we are told there are 3 rings, so we do the display information three 
times. In MACAULAY 2 for a domain, the domain that is the normalization 
is displayed unless we suppress it using a semicolon. When the ring is not a 
domain, a sequence of normal rings, such that the direct sum of these rings 
is the normalization of the reduced ring is displayed. 

Example 2.2. Let A he a union of lines, A = Q[x,y,z]/{{x — y){x — 

z){y - z)). 

11 : A = ClCl[x,y,z] /ideal ((x-y)*(x-z)*(y-z)) 

12 : time integralClosure (A, Variable => V) 

— used 0.47 seconds 

QQ [V , y, z] 

R 

o2 = ( , ) 

2 y - z 

V + V 


o2 : Sequence 

Before proceeding to the Singular example we consider the presentation 
of the first ring. The output of MACAULAY 2 will always give a normal ring, 
but not necessarily a domain. Here it gives 

Q[Vo,y,z] ^ Q[Vo,y,z] Q[Vo,y,z] 
v;^ + Vb ^0 + 1 Vo 

which is normal. In contrast, Macaulay 2 computes the normalization of 
Q[Vb,y, z]/{Vq + Vo) as the isomorphic presentation. 

Using SINGULAR, the output of normal{I) reveals that there are three 
rings. What follows reveals that each of the three rings is isomorphic to 
Q[x,y]. This is the information we gained from Macaulay 2 in an isomor- 
phic format. Since we ran this example in the same session as the previous 
example, we do not need 

> LIB "normal. lib"; 

this time. However, if the Singular session is restarted then this command 
is needed to get the same results. 

> ring R=0, (x,y,z) ,dp; ideal I=(x-y)*(x-z)*(y-z) ; 

> list nor=normal(I) ; 

// 'normal' created a list of 3 ring(s) . 

// nor [3+1] is the delta- invariant in case of choose=wd. 
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// To see the rings, type (if the name of your list is nor): 
show( nor) ; 

// To access the 1-st ring and map (similar for the others) , 

type: def R = nor[l]; setring R; norid; normap; 
// R/norid is the 1-st ring of the normalization and 
// normap the map from the original baser ing to R/norid 

> timer=l; 

> show (nor) ; 

// list, 3 element (s): 
[1] : 

//ring: (0) , (T(l) ,T(2)) , (dp(2) ,C) ; 

// minpoly = 
// objects belonging to this ring: 
// normap [0] ideal, 3 generator (s) 

// norid [0] ideal, 1 generator (s) 

[1] : 

//ring: (0) , (T(l) ,T(2)) , (dp(2) ,C) ; 

// minpoly = 
// objects belonging to this ring: 
// normap [0] ideal, 3 generator (s) 

// norid [0] ideal, 1 generator (s) 

[1] : 

//ring: (0) , (T(l) ,T(2)) , (dp(2) ,C) ; 

// minpoly = 
// objects belonging to this ring: 
// normap [0] ideal, 3 generator (s) 

// norid [0] ideal, 1 generator (s) 

> def Sl=nor[l]; def S2=nor[2]; def S3=nor[3]; 

> setring SI; norid; setring S2; norid; setring S3; norid; 
norid [1]=0 

norid [1]=0 
norid [1]=0 

Since there are three rings whose direct sum forms the normalization of 
A = R/I the output of show(nor); gives the defining rings for all three and 
we see they are the same. Then we ask for the defining ideal and get that it 
is zero, so 

A ^ Q[r(i), r(2)] e Q[r(i), t(2)] e Q[r(i), r(2)] 

which is essentially the same information we got from M ACAULAY 2 . 

Example 2.3. This example was chosen to because this ring has an 

isolated singularity, so it is Ri, but is not S2 in contrast to the first example. 
Thus the normalization of this ring is its S2-ification. This example also 
takes a little longer than the previous examples. Let I be the radical of the 
ideal generated by 

ab^c + bc^d + a^be + c£e + ade^, 
aHc^ + b'^af + a^d^e + ab^e^ + c^de^, 
+ b^ + 0^ + d^ - babcde + 
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We use a few new commands in Macaulay 2 f,his time to make the output 
more presentable for this paper, however these are not needed when a user 
runs Macaulay 2 in emacs since the output is reasonable with wrap-around 
turned off. Even with these new commands we altered the presentation of the 
output a little for both programs to make them more readable in this paper. 

11 : R = QQ[a. .e] ; 

12 : time I = radical (ideal ( 
a*b'"3*c+b*c'"3*d+a'"3*b*e+c*d~3*e+a*d*e'"3 , 
a''2*b*c''2+b''2*c*d''2+a''2*d~2*e+a*b~2*e"2+c~2*d*e~2, 
a"5+b''5+c''5+d''5-5*a*b*c*d*e+e''5) ) ; 

— used 5.77 seconds 

13 : time V = integralClosure (R/I, Variable => X); 

— used 4.23 seconds 

14 : ring ideal V 

o4 = QQ [X ,X ,a,b,c,d,e, Degrees => 

1 {{5},{5},{1},{1},{1},{1},{1}}] 

04 : PolynomialRing 

15 : toString ideal V 

05 = ideal (a"2*b*c"2+b~2*c*d~2+a"2*d"2*e+a*b"2*e"2+c"2*d*e"2, 
a*b''3*c+b*c''3*d+a"3*b*e+c*d"3*e+a*d*e"3, 
a''5+b''5+c''5+d"5-5*a*b*c*d*e+e'"5, 

a*b*c''4-b''4*c*d-X_0*e-a''2*b~2*d*e+a*c''2*d~2*e+b~2*c~2*e~2 
-b*d"2*e"3, 

a*b''2*c''3+X_l*d+a*b*c*d''2*e-a''2*b*e''3-d*e~5, 

a'"3*b'"2*c-b*c"2*d"3-X_l*e-b~5*e-d~5*e+2*a*b*c*d*e'"2, 

a~4*b*c+X_0*d-a*b"4*e-2*b"2*c"2*d*e+a"2*c*d*e"2+b*d"3*e"2, 

X_l*c+b''5*c+a''2*b''3*e-a*b*c~2*d*e-a*d''3*e~2, 

X_0*c-a''2*b''2*c*d-b''2*c''3*e-a''4*d*e+2*b*c*d"2*e"2+a*b*e~4, 

X_l*b-b*c''5+2*a*b''2*c*d*e-c~3*d~2*e+a''3*d*e~2-b*e~5, 

X_0*b+a*b*c"2*d"2-b"3*c"2*e+a*d"4*e-a"2*b*c*e"2+b"2*d"2*e"2 

-c*d*e~4, 
X_l*a-b"3*c"2*d+c*d"2*e"3, 
X_0*a-b*c*d"4+c''4*d*e , 

X_l"2+b"5*c"5+b"4*c"3*d"2*e+b*c"2*d"3*e"4+b"5*e"5+d''5*e~5, 
X_0*X_l+b"3*c"4*d"3-b"2*c~7*e+b~2*c"2*d~5*e-b*c"5*d"2*e"2 
-a*b''2*c*d''3*e''3+b"4*c*d*e~4+a''2*b'"2*d*e"5 
-a*c~2*d"2*e"5-b"2*c"2*e"6+b*d"2*e"7, 
X_0'2+b*c"3*d~6+2*b"5*c*d"3*e+c*d~8*e-b~4*c''4*e~2 

+a~3*c~3*d~2*e~2+2*a~2*b''3*d~3*e''2-5*a*b*c~2*d"4*e~2 

+4*b"3*c"2*d"2*e"3-3*a*d"6*e"3+5*a"2*b*c*d"2*e"4 

-b"2*d"4*e"4-2*b*c"3*d*e"5-a"3*b*e"6+3*c*d"3*e"6-a*d*e"8) 
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The cornmand "ring ideal V" gives the polynomial ring and "toString 
ideal V" gives the defining ideal for the normalization. 

In Macaulay 2 we computed the radical of the three equations and 
then computed the normalization. SINGULAR and Macaulay 2 use differ- 
ent algorithms for the command radical which means that there are examples 
which run well in one program, but not the other. The radical computation 
done above in MACAULAY 2 does not finish in SINGULAR on the Pentium 
III, 600 MHZ machine with 256 MB of RAM, but this is a feature of this 
example. We could just as easily find an example that Singular completes 
and Macaulay 2 does not. In Singular we ran the normalization com- 
putation on the radical obtained from MACAULAY 2. Also this time we do 
not do the command show(nor);, because the output is similar to that of the 
first example. Finally, as we did above, we alter the output a little to make 
it more readable for this paper. 

> ring R = 0, (a,b,c,d,e) ,dp; 

> ideal I = a2bc2+b2cd2+a2d2e+ab2e2+c2de2, 
ab3c+bc3d+a3be+cd3e+ade3 , 
a5+b5+c5+d5-5abcde+e5 , 

a3b2cd-bc2d4+ab2c3e-b5de-d6e+3abcd2e2-a2be4-de6, 

abc5-b4c2d-2a2b2cde+ac3d2e-a4de2+bcd2e3+abe5, 

ab2c4-b5cd-a2b3de+2abc2d2e+ad4e2-a2bce3-cde5, 

b6c+bc6+a2b4e-3ab2c2de+c4d2e-a3cde2-abd3e2+bce5, 

a4b2c-abc2d3-ab5e-b3c2de-ad5e+2a2bcde2+cd2e4; 

list nor = normal (I); 

// 'normal' created a list of 1 ring(s) . 
// nor [1+1] is the delta-invariant in case of choose=wd. 
// To see the rings, type (if the name of your list is nor): 
show( nor) ; 

// To access the 1-st ring and map (similar for the others) , 

type: def R = nor[l] ; setring R; norid; normap; 
// R/norid is the 1-st ring of the normalization and 
// normap the map from the original basering to R/norid 

> timer=l; 

//used time: 3.21 sec 

> def S = nor[l]; setring S; norid; 

norid [1] =T(1) -2*T(2) *T(3) -2+T(2) -2*T(3)*T(4) "2 

+T(1)-2*T(4)"2*T(5)+T(1)*T(2)"2*T(5)"2 

+T(3)'"2*T(4)*T(5)"2 
norid [2] =T(1)*T(2)"3*T(3)+T(2)*T(3)"3*T(4)+T(1)~3*T(2)*T(5) 

+T(3)*T(4)"3*T(5)+T(1)*T(4)*T(5)~3 
norid[3] =T(1) '■5+T(2) "S+TO) -5+T(4) "5 

-5*T(1)*T(2)*T(3)*T(4)*T(5)+T(5)'5 
norid [4] =T(2) *T(3) *T(4) -4-T(3) -4*T(4) *T(5)-T(1) *T(6) 
norid[5]=T(l)*T(2)*T(3) '■2*T(4) '■2-T(2) '■3*T(3) '■2*T(5) 

+T(l)*T(4)"4*T(5)-T(l)"2*T(2)*T(3)*T(5)-2 

+T(2)"2*T(4)"2*T(5)"2-T(3)*T(4)*T(5)~4+T(2)*T(6) 
norid [6] =T(2) -3*T(3) -2*T(4)-T(3)*T(4) -2*T(5) -3-T(l) *T(7) 
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norid[7] =T(1) "2*1(2) "2*1(3) *T(4)+T(2) "2*T(3) "3*T(5) 
+T(1)"4*T(4)*T(5) 

-2*T(2)*T(3)*T(4)"2*T(5)~2-T(l)*T(2)*T(5)-4-T(3)*T(6) 
norid[8] =T(2) *T(3) "5-2*T(l) *T(2) "2*1(3) *T(4) *T(5) 
+T(3)"3*T(4)"2*T(5) 

-T(l) "3*T(4)*T(5) "2+T(2) *T(5) "5-T(2)*T(7) 
norid[9] =T(1) *T(2) *T(3) "4-T(2) "4*T(3) *T(4) 

-T(l) "2*T(2) "2*T(4)*T(5)+T(1)*T(3) "2*T(4) "2*T(5) 

~2*T(3)"2*T(5)"2-T(2)*T(4)"2*T(5)"3-T(5)*T(6) 
)*T(2)"2*T(3)"3+T(1)*T(2)*T(3)*T(4)"2*T(5) 
~2*T(2)*T(5)"3-T(4)*T(5)"5+T(4)*T(7) 
~5*T(3)+T(1)"2*T(2)"3*T(5) 
l*T(2)*T(3)"2*T(4)*T(5) 
l*T(4)"3*T(5)"2+T(3)*T(7) 
~3*T(2)"2*T(3)-T(2)*T(3)"2*T(4)"3-T(2)"5*T(5) 
~5*T(5)+2*T(1)*T(2)*T(3)*T(4)*T(5)"2-T(5)*T(7) 
~4*T(2)*T(3)-T(1)*T(2)"4*T(5) 
(:2)"2*T(3)"2*T(4)*T(5)+T(1)"2*T(3)*T(4)*T(5)"2 
l*T(4)"3*T(5)"2+T(4)*T(6) 

~4*T(3)*T(4) *T(5) "4+2*T(l) "3*T(4) "3*T(5) "4 
:i)"2*T(2)"2*T(4)*T(5)"5 

:i) *T(3) "2*T(4) "2*T(5) "5+2*T(2) *T(4) "2*T(5) "7 
l*T(2)*T(3)*T(4)*T(5)*T(6) 

l*T(3)"2*T(4)"2*T(7)-T(2)"2*T(3)"2*T(5)*T(7) 
:i)"2*T(3)*T(5)"2*T(7)-T(2)*T(4)"2*T(5)"2*T(7) 
l*T(7) 

:i)"2*T(2)*T(4)"4*T(5)"3 

l*T(2)"3*T(4)"2*T(5)"4-T(2)*T(3)"2*T(4)"3*T(5)"4 
(4) "5*T(5) "5+2*T(l) *T(2)*T(3)*T(4) *T(5) "6 

2) "2*T(3)*T(4)"2*T(6)+2*T(1)*T(2)"2*T(5)"2*T(6) 

3) "2*T(4)*T(5) "2*T(6)-2*T(3) "5*T(7) 

4) "5*T(7)+3*T(1)*T(2)*T(3)*T(4)*T(5)*T(7) 

5) "5*T(7)+2*T(7)"2 

i"3*T(2)*T(4)"5*T(5)-T(l)"3*T(3)"3*T(4)"2*T(5)"2 
-2*1(4) "4*T(5) "4+T(2)*T(3) "3*T(4) *T(5) "5 
~2*T(2)"2*T(4)*T(6)+T(2)"2*T(3)"2*T(5)*T(6) 
i:2)*T(4)"2*T(5)"2*T(6)-T(6)"2 

Singular used approximately 1 CPU second less than Macaulay 2 
for this normalization computation. As examples get larger this could he 
significant in obtaining success. 

There are two other operations we want to illustrate. Assume A = R/I 
is an affine domain. These first three examples give the integral closure of A 
in the form S/L where S is of the form i?[yoi • • • ^] and S/L = A. If ^ is a 
domain a user may want the fractions from the quotient field that generate 
the normalization of ^ as a finite algebra over A. This information can be 
obtained using both Macaulay 2 and Singular. 

Example 2.4. For simplicity we use the same hypersurface as Exam- 
vleWli 



+T(2) 
norid[10] =T(1 
-T(l 

norid[ll]=T(2 
-T(l 
-T(l 

norid[12]=T(l 
-T(4 

norid[13] =T(1 
-2*T 
+T(2 

norid[14] =T(2 
+4*T 
-4*T 
+T(1 
+T(1 
-4*T 
+T(6 

norid[15]=2*T 
-T(l 
+3*T 
+3*T 
+2*T 
-3*T 
-2*T 

norid[16]=T(l 
-T(2 
+T(1 
-2*T 
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11 : A = QQ [x,y,z] /ideal(x~6-z"6-y"2*z"4) ; 

12 : time ICf Tactions (A) 

— used 1.77 seconds 

o2 = I x3/z2 x2/z X y z I 

1 5 
o2 : Matrix frac A < frac A 

Thus A[x^ / z'^ , x'^ / z] =A. 

In Singular the method for obtaining the fractions is considerably dif- 
ferent. Several auxiliary functions are used. The library normal. lib is no 
longer used and instead we need reesclos.lib. Singular is very particular 
how you define the ideal defining the affine domain. It must be named ker 
or an error will be issued. The output has the following form. 

fraction; 
[1] : 

xyz8 
[2] : 

yz8 

xyz^ 

Then the fraction we are looking for is [1] divided by [2], that is g- = x. 

The fractions found in this way in Singular are often not reduced. 

We use a new library in this example and since we included the printed 
output that is given when a library is loaded into Singular in Example \2.1\ 
we removed it this time. 

> LIB "reesclos.lib"; 

> ring R=0, (x,y,z) ,lp; ideal ker=x6-z6-y2z4; 

> list L=primeClosure (R) ; closureRingtower(L) ; 

> setring R(6) ; poly f=T(l); closureFrac (L) ; 

> setring R(l) ; fraction; 
[1] : 

xyz8 
[2] : 
yzB 

> setring R(6) ; poly f=T(2); closureFrac (L) ; 

> setring R(l) ; fraction; 
// ** redefining f ** 

[1] : 

y2z8 
[2] : 

yz8 

> setring R(6) ; poly f=T(3); closureFrac (L) ; 

> setring R(l) ; fraction; 
// ** redefining f ** 

[1] : 

yz9 
[2] : 
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yz8 

> setring R(6) ; poly f=T(4); closureFrac (L) ; 

> setring R(l) ; fraction; 
// ** redefining f ** 

[1] : 

x2yz7 
[2] : 

yz8 

> setring R(6) ; poly f=T(5); closureFrac (L) ; 

> setring R(l) ; fraction; 
// ** redefining f ** 

[1] : 

x3yz6 
[2] : 

yz8 

> timer=l; 

//used time: 1.66 sec 

We can use this output from Macaulay 2 and Singular to find a map 
from the variables used by Macaulay 2 to those used by Singular and vice 
versa. Also, given any element of the presentation of A as S/J, Singular 
can find its image in the ring of fractions of A. This is not currently set up 
in Macaulay 2. 

Before computing the integral closure of an ideal, we point out one other 
subtle difference between the two programs. In Macaulay 2 it is possible 
to define multigraded rings and the integralClosure program is designed to 
handle such rings, while in Singular it is not possible to define such rings. 
Multigraded rings play an important role in computing the integral closure 
of an ideal in MACAULAY 2. Since it is not possible to define a multigraded 
ring in SINGULAR the implementation there must be fundamentally different 
from that in MACAULAY 2. 

It is possible, theoretically, to compute the integral closure of an ideal 
using both systems. We say theoretically because the computation is often 
much too complex to complete, either due to memory or time (mostly mem- 
ory). The approach, in both systems, is the classical one, that is compute 
the Rees algebra first, find the normalization of the Rees algebra and then 
find the degree one piece of that graded algebra. 

In Macaulay 2 a function to compute I is not yet implemented in 
the main distribution, but we include one in the appendix which is being 
submitted to Macaulay 2 for inclusion. To run this program, a program to 
compute the Rees algebra that preserves the natural multigrading is needed. 
We also include the code for this in the appendix. There are several such 
programs for MACAULAY 2 in circulation. Besides the one in the appendix, 
one can be found in jl4j and we have received yet another via personal 
communication from David Eisenbud. A package of such programs is being 
put together for inclusion with MACAULAY 2. The one we include here is 
the only one that incorporates the natural multigrading, but is otherwise 
fundamentally the same as that in jl4j . which could be easily altered to use 
the grading. The program communicated by David Eisenbud is more general 
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and is based on his paper with Huneke and Ubich Both computing 
blowups and integral closure of ideals are functions in Singular if the Rees 
library has been loaded. 

We include two complete examples and discuss briefly two others. The 
first example is one that can be checked by hand using standard techniques. 
The second example is considerably more complicated and the remaining 
ideals illustrate how easily this process fails to complete. All the examples 
are small in the sense that they use at most three generators, use at most 
three variables and have degree at most seven. 

Example 2.5. Let R = Q[x,y] and I = {x^ ,x^y^ ,y^). Using the fact 
that the integral closure of a monomial ideal corresponds to the integer lattice 
points in the convex hull of the exponent vectors of the ideal, this example 
is easily computed by hand. The first two commands load files named Ideal- 
Normal. m2 and blowup. m2 which are the programs given in the appendix. 

11 : load "IdealNormal .ni2" ; load "blowup .ni2" ; 
— loaded IdealNormal .m2 

— loaded blowup. m2 

12 : R = QQ[x,y]; I = ideal(x"2,x*y"4,y"5) ; 

03 : Ideal of R 

14 : time ideallC(I) 

— used 0.98 seconds 

2 3 5 

04 = ideal (x , x*y , y ) 

o4 : Ideal of R 

Computing the integral closure of an ideal in Singular uses the library 
reesclos.lib which we used in Example \2.4\ The library must be loaded if it 
has not already been. 

> ring R = 0,(x,y),dp; ideal I = x2,xy4,y5; 

> list J = normall(I); J; 
[1] : 

_[l]=x2 
_[2]=y5 
_ [3] =-xy3 

The output following [1]: is the generators of I. 

Example 2.6. Let R = Q[x,y, z] and I = {y^ + x'^z, —x^ + y^z'^). 

15 : R = QQ[x,y,z, MonomialSize => 16]; 

16 : I = ideal (y"6+x"2*z, -x"6+y~4*z~2) ; time ideallC(I) 



o6 : Ideal of R 

— used 3.6 seconds 
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6 2 6 4254 23 
o7 = ideal (y +xz,x -yz,xy + x*y z ) 

o7 : Ideal of R 

> ring R=0, (x,y,z) ,dp; ideal I=y6+x2z,-x6+y4z2; 

> list J=normalI(I) ; J; 
[1] : 

_[l]=y6+x2z 
_ [2] =-x6+y4z2 
_ [3] =x5y4+xy2z3 

These examples show that both Macaulay 2 and Singular can com- 
pute the integral closure of a 2 generated ideal in 3 variables, but an ideal 
with three generators will often cause problems. For example, if i? = 
Q[x,y,z] and I = {y^ + x^z,-x^ + y^z^.x^y - z^y) then Macaulay 2 
will give a monomial overflow error and Singular runs for about 7.1 hours 
and uses up the memory available on this machine. But, with two variables 
and two generators we can find a troublesome ideal in degree as low as seven. 
If i? = Q[x,y] and / = (y'^ + x^, —y^ + xy'^) then Macaulay 2 finishes in 
905.29 CPU seconds and uses about 66% of the memory. The ideal we get 
is 

7 = {x^y^ + x'^, x^y^ - x^ x^y + 2,9x'^y'^ - 38x^, x^ - x^ y^ - xy^, xy^ + x^). 

We ran this same example in Singular and after 80,836.49 CPU seconds 
the machine quit the computation after using all of the memory. Both 
programs struggle with various computations of ideal integral closure due to 
computing the integral closure of the blow-up ring which has more variables 
than the original ring. Thus, while it is possible to compute the integral 
closure of some ideals using this method, for most cases it is impractical and 
the myriad of papers dealing with special cases must be consulted. 

3. Appendix 

The program ideallC, below, computes the integral closure of an ideal 
in Macaulay 2. This is the program used to compute the closures in 
Examples 12. 5( 12.61 Except for ICfractionsLong all of the functions used 
in ideallC are self-explanatory, are commented, or are clearly explained in 
the Macaulay 2 documentation. ICfractionsLong is similar to the function 
ICfractions used in Example 12.41 which computes the normalization of A and 
then returns a minimal generating set for A as an algebra over A. During the 
normalization computation fractions are computed that may be extraneous 
and while ICfractions only returns a minimal set, ICfractionsLong keeps 
track of all the fractions generated by the computation. To ensure that the 
entire degree one component of the integral closure of the Rees algebra is 
computed, ICfractionsLong is used rather than ICfractions. 
ideallC = methodO 

ideallC (Ideal) := Ideal => (I) -> ( 
R := ring I; 
nl := numgens R; 

Jl := blowup (I) ; — defining ideal of blow-up 
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R2 := ring Jl; 

ii2 := numgens R2; 

S := R2/J1; ~ the blow-up 

Sfracl := first entries ICf ractionsLong S; 

— The slow part is ICf ractionsLong(S) , as expected. 
Sfrac2 := apply (#Sfracl-n2, i-> Sfracl#i) ; 

toLift := select (Sfrac2 , i-> (degree i)#l == 1); 

— toLift is the set of elements in the normalization 

— of the Rees algebra that generate the degree one 

— component . 

Varlmage := flatten append (gens R, first entries gens I); 

— We use Varlmage to define the map from R2 -> R needed 

— to lift the fractions back to R. 
LiftMap := map(R, R2, Varlmage) ; 
NewNums := apply (toLift , i-> LiftMap ( 

substitute (numerator i, R2))); 
NewDenoms := apply(toLif t , i-> LiftMap ( 

substitute (denominator i, R2))); 
NewGens := apply (#toLift , i-> substitute ( 

(NewNums#i)/ (NewDenoms#i) ,R)) ; 
ideal mingens (I + ideal (NewGens)) 
) 

A program for computing the Rees algebra R[tl] in Macaulay 2 which 
is used in the program IdeallC give above. 

blowup = method(Options => {VarName => Y}) 
blowup (Ideal) := o-> (J)-> ( 

— Input: J is any ideal. 

— Output: The result is the defining ideal of the 

blowup algebra R[tJ] where R is the ring 
of J. 

— METHOD: We construct a polynomial ring 

R[t,y_l, . . ,y_n] where n is the nimber of 
generators of J. We then construct the ideal 
(Y_l-tJ_l, . . . ,Y_n-tJ_n) . Then R[tJ] is 
isomorphic to the polynomial ring 
R[Y_1, . . . ,Y_n] mod the out put of this 
algorithm. 

R := ring(J) ; 

nl := numgens(J); 

n2 : = numgens (R) ; 

Degsl := flatten degrees source gens J; 

— Degsl are the degrees of the generators of J. 
Degsla := apply (nl , i->l+Degsl_i) ; 

Degs2 := flatten ((monoid R) . Options .Degrees) ; 

— Degs2 is the degrees of the variables of R. 
Degs := join({l}, Degs2, Degsla) ; 

S := coeff icientRing(R) [tjgens R, 

(o.VarName)_(l) . . (o. VarName) _(nl) , 
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MonomialOrder => Eliminate 1 , Degrees => Degs] ; 
Jl := substitute ( J, S) ; — Put J in the bigger ring S. 
M := inatrix{apply(nl,i-> S_ (i+n2+l) -S_0*Jl_i) } ; 

— This is the matrix (Y_l-t J_l , . . . , Y_n-t J_n) 
K := gb(M); 

— Next eliminate t. 

L := ideal (selectlnSubringCl , gens (K) )) ; 

— Set the Multidegrees . 

Degslb := apply(nl,i->join({Degsl_i},{l})) ; 

Degs2b := apply (n2 , i-> join({Degs2_i> , {0}) ) ; 

Degsb := join(Degs2b, Degslb) ; 

Sb := coef f icientRing(R) [gens R, 

(o . VarName) _1 . . (o . VarName) _ (nl) , 
MonomialOrder => Product0rder{n2,nl}, 
Degrees => Degsb] ; 

trim(substitute(L,Sb) ) 

) 
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